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We have measured the bias voltage and field dependence of eigenmode frequencies in a magnetic 
tunnel junction with MgO barrier. We show that both free layer (FL) and reference layer (RL) 
modes are excited, and that a cross-over between these modes is observed by varying external field 
and bias voltage. The bias voltage dependence of the FL and RL modes are shown to be dramatically 
different. The bias dependence of the FL modes is linear in bias voltage, whereas that of the RL 
mode is strongly quadratic. Using modeling and micromagnetic simulations, we show that the linear 
bias dependence of FL frequencies is primarily due to a linear dependence of the perpendicular spin 
torque on bias voltage, whereas the quadratic dependence of the RL on bias voltage is dominated 
by the reduction of exchange bias due to Joule heating, and is not attributable to a quadratic 
dependence of the perpendicular spin torque on bias voltage. 
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I. INTRODUCTION 

A current flowing in a magnetic structure will be spin- 
polarized because of spin-asymmetric scattering in the 
two spin channels (we shall assume that the spin relax- 
ation length is longer than any relevant system dimen- 
sion). In non-collinear systems or magnetic heterostruc- 
tures in which the magnetization in different layers is 
aligned in different directions, the spin polarized conduc- 
tion electrons can exert a torque on the magnetization or- 
der parameters. This spin transfer torque couples the di- 
rect current with magnetization dynamics and allows for 
the manipulation of magnetization using spin-polarized 
currents.^* 2 - In a magnetic multilayer system, the spin 
transfer torque can be used to pump energy into the mag- 
netization dynamics so as to precisely offset dissipative 
losses, leading to self-sustaining magnetization oscilla- 
tions. These so-called spin torque oscillators (STOs) can 
potentially be used for microwave signal generation^ 
modulation;^ and detection^ Recently, STOs based on 
CoFeB/MgO/CoFeB magnetic tunnel junctions (MTJs) 
have attracted considerable interest because of their rel- 
atively large microwave power In addition to the in- 
plane spin torque predicted by Slonczewski and Berger^^ 
the perpendicular spin torque first predicted by Zhang, 
Levy, and Ferfc^i is appreciable in MTJ o 12 i 13 and much 
larger than in metallic systems*^— As the perpendicu- 
lar spin torque couples to the magnetization dynamics, 
it is of significant interest to determine this component 
for fundamental understanding as well as for future ap- 
plication*^ 

The perpendicular spin torque in an MTJ depends on 



the current, and therefore on the bias voltage applied 
across the MTJ. Several methods have been proposed 
and used to determine the magnitude and bias depen- 
dence of the perpendicular spin torquei 8 ' 12 i 13 i 18 ~ ^ The 
reported results vary both in their qualitative and quan- 
titative assessment of the bias dependence of the per- 
pendicular spin torque. Some studies show a linear ) 18 i 23 
others show a quadratic , 12 d 9 i 20 i 24 or a combination of lin- 
ear and quadratic^! behavior of the perpendicular spin 
torque with bias voltage. Early theoretical works pre- 
dicted a quadratic dependence of spin torque on bias 
voltage . 25 ! 26 A recent theory however predict a linear- 
contribution to perpendicular spin torque for asymmet- 
ric MTJs, in which the two magnetic layers on either 
side of the tunnel barrier are not identical, which was 
the case for Ref. [U However, in a recent study of some 
400 MTJ devices combined with micromagnetic simula- 
tions^ 3 - it was shown that the perpendicular torque varies 
linearly with bias voltage (for a certain range of bias volt- 
ages) even for nominally symmetric MTJs, and that any 
quadratic term is negligible. 

We will here present results for the bias and field de- 
pendence of eigenmode frequencies for MTJs with MgO 
tunneling barriers and nominally symmetric interfaces on 
each side of the tunneling barrier. In reality, there will 
be some asymmetry introduced by deposition kinetics 
and thermal anneal, but even large, deliberately intro- 
duced asymmetry have been shown to have a relatively 
small effect on the perpendicular spin torque^ and we 
do not expect that asymmetries play any significant role 
in the behavior of the MTJs studied here. A cartoon 
of the cross-section of the devices is shown in the in- 
set of Fig. [U An antifcrromagnctic IrMn layer provides 
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Figure 1. (Color online) Magnetoresistance loop of the MTJ 
nanopillar showing TMR of 70%. The top inset shows a mag- 
nified view of the low field region whereas the bottom inset 
shows the layer stack of the MTJ. 



an exchange bias on the CoFe pinned layer (PL) which 
tends to keep its magnetization direction in a fixed di- 
rection. The PL is strongly coupled antifcrromagneti- 
cally through a 0.81 nm thick Ru layer to the composite 
CoFe/CoFeB fixed layer, or reference layer (RL). This 
CoFe/Ru/CoFe/CoFeB structure is referred to as a syn- 
thetic antiferromagnet (SAF) as its net magnetic mo- 
ment is close to zero. Above the MgO tunnel barrier is 
the CoFeB free layer (FL) which can easily be rotated 
by an external field. We will show that FL as well as 
RL magnetization cigenmodes can be observed, depend- 
ing on the strength of the applied magnetic field. While 
the observation of magnetization modes in the GHz range 
has been already reported for similar devices^ - — the de- 
tailed bias dependence of the FL and RL modes and the 
role of perpendicular spin torque on these modes have not 
been discussed. For the first time, we show that the bias 
dependence of the FL and RL modes are dramatically dif- 
ferent. The bias dependence of the FL modes is approxi- 
mately linear in bias voltage, whereas the RL modes show 
a strong parabolic dependence with bias voltage. Using 
modeling and micromagnctic simulations, we show that 
this linear bias dependence of the FL frequencies is pri- 
marily due to a linear dependence of perpendicular spin 
torque, whereas the parabolic decrease of frequency of 
the RL with bias voltage is dominated by the reduction 
of exchange bias due to Joule heating. 



II. EXPERIMENT 

The MTJ samples were fabricated using state-of- 
the-art industrial ultra-high-vacuum sputter deposi- 
tion.— £L21 The deposition technique typically produces 
layers with a root mean square roughness of about 0.1 nm 
after thermal anneal. We will here predominantly discuss 
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Figure 2. (Color online) (a) Mode frequencies versus easy 
axis external field measured at low bias current of 1 mA. The 
dashed lines define the boundary of the three regions, I, 77 and 
777 as described in the text, (b) Example spectra measured 
at 77 cxt = 400 Oe and low bias current of 7b = 1 mA. 



room temperature results from a circular device with an 
approximate diameter of 240 nm. However, similar re- 
sults are also obtained in a number of other devices of 
same dimensions, as well as several devices with diameter 
of 180 nm. The magnetization densities of the FL and 
RL are approximately 1000 emu/cm 3 , and that of the 
pinned layer approximately 1200 emu/cm 3 . 23 i 33 ' 34 The 
direction of the exchange bias in all devices is defined 
as the x-direction. The measurement set-up is similar 
to that described in Ref. The signal generated from 
the STO was amplified using a broadband 0.1-26 GHz, 
+45 dB microwave amplifier, and detected by a spectrum 
analyzer. The dc bias current is fed to the device by a 
current source through a 0-26 GHz bias-tee connected in 
parallel with the transmission line. We use the conven- 
tion that a positive current flows from the FL to the RL, 
in which case the electrons flow from the RL to the FL. 

In Fig. [TJ a magnetoresistance curve is shown where 
positive field is along the x-direction which is also the 
direction of the exchange bias on the PL. The measured 
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Figure 3. (Color online) Bias dependence of mode FLAP1 
at f/ext = 400 Oe and a linear fit (red solid line). The inset 
shows power spectral density (PSD) of this mode as function 
of bias voltage. 
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Figure 4. (Color online) Measured bias dependence of the 
mode RLAP1 at H ex t=1000 Oe (solid squares) and calcula- 
tions using Eqs. © and ([SJ for bj = bi = (blue solid line), 
and bj — b\Vh with bi — 37 Oe/V (red dashed line). The inset 
shows power (PSD) of this mode as function of bias voltage. 



III. RESULTS 



tunneling magnetoresistance (TMR) is about 70% and 
the resistance area (RA) product in the parallel state is 
about 1.5 fi (fim) 2 . We found the expected strong bias 
dependence of TMR and decrease of resistance with tem- 
perature, both of which confirm the integrity of tunneling 
barrier and absence of shorts * 3 2 The temperature de- 
pendence of the TMR in our devices is measured to be 
-0.13%/K. The bias range was limited to \V h \ < 0.4 V to 
avoid any risk of damage to the MTJs due to prolonged 
exposure to bias voltages. During the measurement, we 
continuously measure resistance of the device to make 
sure that the device is not shorted. After application of 
high voltages (|Vb| > 0.6 V), some samples exhibited a 
strong reduction of TMR similar to Ref ITol Such low 
TMR samples are not investigated in this work. We ob- 
serve a loop shift of about 100 Oe of the R — H curve (as 
shown in the inset) which is a result of RL-FL ferromag- 
netic interlayer exchange coupling (IEC) due to, e.g., bar- 
rier roughness (Neel orange-peel coupling) or electronic 
processes such as coherent scattering processes that also 
give rise to an RKKY-like coupling. The coercive field 
of the FL is approximately 20 Oe. At higher fields, the 
resistance decreases due to rotation of the SAF. The two 
SAF layers (PL and RL) starts to rotate when the net 
field on the PL exceeds the exchange bias field. We ob- 
serve a small hysteresis in this regime indicating that the 
rotation of the two layers is not fully coherentj 43 ' 44 From 
the rotation of the SAF, we define the exchange bias field 
.ffeb using the peak in the first derivative of resistance ver- 
sus field on the positive field side. The magnetoresistance 
loop of Fig. [T] was measured at a low bias I\, = 0.1 mA, 
from which we found i? c b = 2200 Oe. 



In Fig. [Ha), we show the low-bias behavior of fre- 
quency versus external field, H cxt . A bias current of 
lb = 1 mA was used for this measurement (which corre- 
sponds to about 42 mV in the P state and 73 mV in AP 
state) . The frequency of the modes was determined from 
the measured spectra. An example spectra is shown in 
Fig.OHb) for H ext = 400 Oe and J b = 1 mA. Based on the 
low bias measurement of Fig. [Ha), we distinguish three 
regions in which three different behaviors of the eigen- 
modes are observed. The first region (J) corresponds to 
positive field, 100 < H cxt < 800 Oe in which the FL and 
RL magnetizations are in an antiparallel (AP) configura- 
tion. In this region the frequencies of the modes increase 
with external field. As we will demonstrate later, these 
modes correspond to free layer excitations and we will 
label them FLAP. The second region (77) corresponds to 
higher positive fields, H oxt > 800 Oe, where additional 
lower-frequency modes appear with frequencies that de- 
crease with -f/ext ■ We will identify these new modes as RL 
modes, and we will label them RLAP. Finally, the third 
region (///) corresponds to the case of H ext < 100 Oe 
with the magnetization of the FL and RL in parallel (P) 
configuration. The frequency of the single mode visible in 
this region increases with external field strength. We will 
identify this mode too as a FL mode. To distinguish from 
the FL modes in the AP configuration we will label it as 
FLP1. To study the bias dependence we identify three 
modes in these three regions labeled as FLAP1, RLAP1, 
and FLP1 as shown in Fig. dta), which arc strongest 
(with highest power) modes in the respective regions. 

The bias dependencies of the frequency and power of 
the FLAP1, RLAP1, and FLP1 modes are shown in 
Fig. [31 Fig. 0] and Fig. [SJ respectively. The field val- 
ues were chosen such that apparent mode-crossings are 
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Figure 5. (Color online) Bias dependence of mode FLP1 at 
Hext = —400 Oe and a parabolic fit (red solid line). The inset 
shows power (PSD) of this mode as function of bias voltage. 



age, but shows the opposite bias asymmetry, i.e. it is 
higher for electrons flowing from FL to RL, indicating 
that the excitation is a RL mode. As will be discussed in 
the next section, this is also confirmed by micromagnetic 
simulations and modeling. 

The frequency of the FLP1 mode in Fig. [5] also shows 
a parabolic dependence on bias voltage, which is again 
in contrast to that of the FLAP1 mode. However, the 
parabolic dependence is much weaker than that of the 
RLAP1 mode with a curvature that is only about 30% of 
the RLAP1 mode's curvature. The power of this mode 
also increases with bias voltage. However, the maxi- 
mum power is by about one order of magnitude lower 
compared to the RLAP1 mode and about two orders of 
magnitude lower compared to that of the FLAP1 mode. 
Interestingly, the power has a polarity-dependent asym- 
metry similar to that of the RLAP1 mode. 

IV. DISCUSSION 



avoided and there is minimum mixing of modes from dif- 
ferent branches. Firstly, the FLAP1 mode shows a linear 
behavior of frequency with bias voltage for |Vb| < 0.4 V. 
The power spectral density shows a strong polarity- 
dependent emission power and increases almost exponen- 
tially with bias voltage. Since the structure is in the AP 
state at £f ex t = 400 Oe, spin torque tends to destabilize 
the FL magnetization for positive polarity. Hence higher 
power for positive polarity is consistent with spin torque 
excitation of FL layer. 

The frequency of the RLAP1 mode on the other hand 
shows a strong parabolic behavior as a function of bias 
voltage (Fig.HJ), and a negligible linear contribution. The 
bias dependence of the frequency is hence symmetric in 
bias voltage in contrast to that of the FLAP1 mode, 
which is linear in bias voltage and hence antisymmet- 
ric with Vb (apart from a trivial constant). The power of 
the RLAP1 also increases exponentially with bias volt- 

I 



A. Low-bias mode frequencies 

We start by analyzing the low-bias behavior of the 
FLAP1, FLP1 and RLAP1 modes shown in Fig. Ufa). 
We take a coordinate system with x along the positive 
field direction (along the direction of the exchange-bias 
field on the pinned layer), y perpendicular to the field 
and in the plane of the magnetic layers, and z perpen- 
dicular to the plane of the layers. A simple Kittel-like 
model can be obtained for the lowest bulk-like FLAP1 
mode (we are for now ignoring edge modes) in the fol- 
lowing way. We represent the magnetization density as 
a macro-spin vector Ms(t) = Msm(t), where Ms is the 
saturation magnetization density of the FL and m(t) a 
unit vector along the magnetization direction. If we as- 
sume that stray fields from the SAF can be ignored, the 
effective field acting on the FL in the AP configuration 
with the RL magnetization along — x is 



H cff = (H cxt - H mc ) x ~ M s (N x m x x + N y m v y) - bjy - H d m z z. 



(1) 



Here -ffiEc is the effective field due to the RL-FL IEC, 
N x and N y are the xx and yy components of the demag- 
netizing tensor (assuming an ellipsoidal shape for which 
the demagnetizing tensor is diagonal), H e i = AttMs is 
the out-of-plane demagnetizing field, and bj is the ef- 
fective field due to perpendicular spin torque. We have 

I 



I 

here included bj for completeness, but in the low-bias 
limit we will later set bj = 0. By linearizing the torque 
equation dMs/dt — ~\j e \Ms x -Hcff in small excur- 
sions about the equilibrium magnetization configuration, 
m(t) = x+Sm(t), with \5m\ <C 1 and perpendicular to x, 
we obtain for the zero-bias resonant frequency wflapi.o 



•>flapi,o = ll [(Ny - N x ) M s + # ext - Hjec - bj] [H d - N X M S + H cxt - H mC - bj] , 



(2) 



I 

with 7 e is the electron gyromagnetic factor, \y e \ « prctation: apart from the demagnetizing fields, the main 
2.8 GHz/kOc. Equation [2] lends itself to a simple inter- 
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Figure 6. (Color online) Fit to the FLAP1, RLAP1, and FLP1 
modes of Fig. [2|a) using Eqs. (|2| - Q 
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Figure 7. Experimental values (filled circles) for the modes 
FLAP1 (black), RLAPf (blue), and FLP1 (black), together 
with micromagnetically calculated frequencies for the lowest- 
lying bulk modes (squares, triangles, diamonds, and circles). 
The agreement is satisfactory, and the micromagnetic results 
show a cross-over to a branch with decreasing frequency as 
a function of magnetic field, similar to the experimentally 
measured cross-over from the FLAP1 to the RLAP1 modes. 



stiffness and restoring torque is provided by the exter- 
nal field, while the ferromagnetic IEC softens the mode, 
since the RL magnetization is in the —x direction. In 
the P configuration, we similarly obtain for the zero-bias 
FLP1 mode (assuming the RL is stationary) 



<4lpi,o = ll i( N v ~ N*) M s + #cxt + Hiec + bj] [H d - N X M S + H cxt + H 1EC + bj] 



(3) 



Finally, again under the assumption that we can ignore stray fields from the SAF, we obtain for the zero-bias 
RLAP1 mode 



Wrlapi.o = 7e [N X M S - H cxt + H ch - H 1EC + bj] [H d + N X M S - H eXt + H ch - H mC + bj 



(4) 



For this mode, the external field softens the mode as 
it is antiparallel to the equilibrium direction of the RL 
magnetization, while the IEC stiffens the mode. It is 
now clear that Eqs. |[2J - Q at least qualitatively can 
explain the behavior of the FLAP1, RLAP1, and FLP1 
modes: The cross-over from FLAP1 to RLAP1 modes 
occurs when the external field is large enough (for a 
given IEC) that it drives the RLAP1 frequency below 
that of the FLAP1 mode. Figure [6] shows a fit to the 
data in Fig. [5Ja) (positive bias) using Eqs. ©-(HI)- For 
the FLAP1 and FLP1 modes the same parameters were 
used, with (N y - N X )M S = 100 Oe, Hiec = 100 Oe, and 
H d ~ N X M S = 7 kOc, and bj = 0. To fit RLAP1 mode, 



we used Eq. H with N X M S = -15 Oe, H ch = 2.2 kOe, 
and H d + N X M S = 4.0 kOe^^^ The values for the 
FLAP1 and FLP1 modes are quite reasonable, especially 
because the dimensions of the device have not been de- 
termined accurately, but only the target dimensions from 
the lithographic process are known. Also, the magnetiza- 
tion densities for the layers are estimated based on values 
from similar devices and process conditions. The values 
used for the RLAP1 fit are also quite reasonable in view 
of these facts. This simple analysis supports our identi- 
fication of the modes as FL modes (FLAP1 and FLP1), 
and a RL mode (RLAP1), and therefore also the cross- 
over of the lowest-lying mode from FLAP1 to RLAP1 at 
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Figure 8. (Color online) Eigenmode distributions for the two 
lowest eigenmodes at H eyi t = 400 Oe (left panel) and H BXt — 
1000 Oe (right panel). In the figure, the distribution of the 
real part of the x-component of the eigenmode is plotted for 
the FL at 400 Oe and the RL at 1000 Oe. At each of these 
field values, the exhibited distribution of the eigenmodes was 
dominant, with distributions in the RL (FL) negligible at 400 
Oe (1000 Oe). The figure clearly shows a cross-over from a 
FL bulk-like mode, to a bulk-like RL mode. 



H cxt « 800 Oe. 

To further support this analysis, we also performed mi- 
cromagnetic modeling of a 240 x 240 nm 2 device with 
layer thickness PL=2 nm, RL=3 nm, and FL=3 nm, 
and magnetization densities as described above; the IEC 
coupling was taken as 0.035 erg/cm 2 , corresponding to 
Hmc ~ H7 Oe, the pinning energy to the AFM as -1.2 
erg/cm 2 ; the intralayer exchange coupling was taken as 
1.3 x 10 -6 erg/cm, and the dimensionless damping con- 
stant a = 0.1i 23 i 33 ' 34 The system was divided into 5 nm 
x 5 nm cells in the xy-plane and the magnetization con- 
figurations obtained by integrating the time-dependent 
Landau-Lifshitz-Gilbcrt (LLG) equations with demagne- 
tizing fields calculated using fast Fourier transforms in 
the xy-plane and direct sums in the z-direction. The 
system was initialized with random magnetization direc- 
tion in each cell, then subjected to an external field of 
4 kOe along the x-direction and the magnetization equi- 
librated. The field was then changed to 1.4 kOe and 
subsequently reduced in steps of 200 Oe to -1 kOe, and 
the magnetization equilibrated at each field value. At 
each new field value the magnetization was given a small 
random perturbation from its previous equilibrated con- 
figuration. At each equilibrated magnetization config- 
uration, the LLG equations were linearized about the 
equilibrium and the resulting homogeneous linear equa- 
tions solved for the (complex) magnetization eigenmodes. 
Figure [7] shows the experimentally measured frequencies 
for FLAP1, FLP1, and RLAP1 modes together with the 
real part of the calculated frequencies for the lowest- lying 
bulk modes (some modes that were clearly edge modes 
were excluded from the plot). The qualitative behavior 
of the calculated frequencies agrees with that of the mea- 
sured ones, and the quantitative agreement is reasonable. 
The main point here is that the calculated frequencies 
for these branches all decrease with increasing field from 



about 800 Oe. Detailed examination of the distribution 
of magnetization of these branches confirm that the ex- 
cited magnetization motion shifts from the FL at 400 
Oe to the RL at 1000 Oe (Fig. [5). This cross-over from 
FL to RL mode is rather insensitive to the parameters 
used and is also observed for different geometries, e.g., 
300 x 600 nm 2 , 150 x 600 nm 2 or 150 x 300 nm 2 . We also 
note that the FLP1 branch of these three modes (nega- 
tive external field) tends to be distributed spatially near 
an edge on the x-axis, although the precise distributions 
depend on the particular parameters used. This is not 
inconsistent with the experimental observation that the 
power of the FLP1 mode is much lower than that of the 
FLAP1 and RLAP1 modes. 



B. Bias- voltage dependence of the mode 
frequencies 

We now turn to the bias voltage dependence of the fre- 
quencies of the FLAP1, RLAP1, and FLP1 modes. We 
can safely ignore any bias dependence of magnetization 
density, IEC, and intralayer exchange, and will focus on 
the bias dependence of the perpendicular spin torque ef- 
fective field bj and the exchange bias field rJ e t>- Also, 
we note that at these field values the FL and RL re- 
main in P or AP configuration for relatively small bias 
voltages under consideration, so we take the resistance 
for each configuration to be constant and independent of 
bias voltage. Under the assumption that the change in 
bj and .Heb with Vb are small compared to their zero-bias 
values, we can expand bj and H cb in Taylor series about 
Vb = 0. We will ignore any observable direct spin torque 
effect on the exchange bias 4^ It is well established that 
the bias dependence of H e \> decreases monotonically with 
temperature (and has to vanish at the Neel temperature 
of the AFM) i 46 i 47 For our samples, we measure the tem- 
perature dependence to be is -1.65 Oe/K, which is consis- 
tent with literature for a similar dcvice4£ Therefore, the 
main bias voltage dependence of H e b comes from Joule 
heating, and H b must then be an even function of Vb, 

H ch (V h ) = H° h -aV*, (5) 

where H® h is the zero-bias value of H e \, and a is a con- 
stant of dimension OcV~ 2 . In contrast, we will for now 
expand the perpendicular spin torque in both linear and 
quadratic terms, 

bj(V h ) = b° J + b 1 V b + b 2 v£, (6) 

where b°j is the zero-bias value of bj, and the constants 
&i and &2 have the dimensions of OeV -1 and OeV~ 2 , 
respectively. We will assume that the constants are the 
same for both AP and P magnetization configurations^. 
Inserting Eqs. ([3]) and ((BJ) into Eqs. ©-Q and expanding 
the square roots to second order in Vj, we obtain 
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Figure 9. (Color online) Bias dependence of H E h and a 
parabolic fit (red solid line). The inset shows a linear fit of 
H e b to applied power. 
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We first discuss the bias dependence of the FLAP1 
mode. If we compare Eq. ([7]) with Fig. [31 it is clear that 
the linear dependence of the frequency of this mode on Vb 
implies that &2 ~ so that the perpendicular spin torque 
effective field is linear in Vb- This is in agreement with 
the results of Petit et al.^ and Hcinoncn et al££. Using 
the values for (Ny - Nx)M s = 100 Oc, H mo = 100 Oc, 
ff ext = 400 Oe, and H d - N X M S = 7 kOe from above, 
we can estimate the constant b\ as b\ ~ 37 Oe/V. 
With RA = 1.5 ilptm 2 we then obtain a perpendicular 
spin torque field that is about 0.55 x 10~ 6 Oc(cm) 2 /A. 
This is the same order of magnitude as the value of 
10 -6 Oc(cm) 2 /A obtained by Petit et alM on A1 2 3 
MTJs with a TMR of 20%, but an order of magnitude 
smaller than the value of about 11 x 10~ 6 Oe(cm) 2 /A 
obtained by Heinonen et alM on MgO MTJs with RA « 
0.5 fi/Um 2 and a TMR of about 85%. It is not clear why 
there is a relatively large spread in these observed values. 
Tight-binding^ and s-d models^ suggest that the spin 
torque is sensitive to the particular parameters used in 
the model (e.g.., band offsets, s-d coupling) as well as to 
how much of the transverse spin density is absorbed in 
the FL. We speculate that the observed large variation in 
the perpendicular spin torque is a reflection of detailed 
differences in the MTJs used here and in Refs. [HI andl23l. 

Next, we turn to the bias dependence of the RLAP1 
mode, Eq. It is clear from this equation that a re- 
duction in H c h due to Joule heating can have a dominant 



effect on the bias dependence of this mode's frequency. In 
order to clarify this effect, we measured H e h as a function 
of Vb (Fig. |9j), from which we obtain H eh = H° h - aV^, 
with a = 1880 Oc/V 2 . We insert this with the same val- 
ues as those used in Fig. [5] together with b\ =37 Oe/V 
in Eq. © to get the result shown in Fig. SJ The scale of 
the exchange bias field completely dominates the voltage 
dependence of this mode's frequency compared to the ef- 
fect of the spin torque. In Fig. @] is also exhibited the 
frequencies obtained by setting b\ = 0. The result is al- 
most indistinguishable from that obtained by including 
the perpendicular spin torque field. We therefore con- 
clude that the parabolic dependence of this mode's fre- 
quency on bias voltage is entirely explained by the weak- 
ening of the exchange bias field due to Joule heating, and 
is not due to spin torque. 

In the analysis of the FLAP1 and RLAP1 modes, we 
have ignored the fact that the resistance in the AP state 
decreases linearly with increasing magnitude of bias volt- 
age^ In the present case, the AP resistance decreases 
at a rate of about 30 f2/V. For the FLAP1 mode, we 
could use current instead of bias voltage as independent 
variable and consider the change in frequency from its 
zero-bias value as a function of current. This gives a 
slope of about -17 GHz/ A, which is equal to the slope 
as a function of bias voltage, -0.23 GHz/V, multiplied 
by the zero-bias resistance, 74 f2, to within experimental 
uncertainties. Similarly, we can correct for the change 
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in resistance in the analysis of the RLAP1 mode. How- 
ever, such a correction does not alter our findings in any 
qualitative or quantitatively significant way. 

Finally, we turn to the bias dependence of the FLP1 
mode, which shares many of the characteristics of the 
RLAP1 mode. As we noted earlier, the bias depen- 
dence of this mode's frequency is parabolic, similar to 
the RLAP1 mode, however with a weaker curvature. The 
power of the FLP1 mode also shows the same bias volt- 
age asymmetry as the RLAP1 mode, i.e. it is stronger 
for electrons flowing from the FL to the RL, which is 
somewhat surprising for an FL mode. It also exhibits a 
much lower power than either of the RLAP1 or FLAP1 
modes. Micromagnetic modeling suggests that this mode 
is a FL mode, but the precise spatial distribution (edge 
vs. bulk) depends rather sensitively on the specific pa- 
rameters used in the modeling. To lowest order, the FL 
modes are not sensitive to the exchange bias strength, so 
one possibility is that the parabolic bias dependence ob- 
served for this mode is indeed due to perpendicular spin 
torque. A fit to the experimental data in Fig. [S] gives 
/ - fo ~ -2.8 GHz/V 2 , which inserted in Eq. © gives 
fe 2 ~ —3 x 10 2 Oe/V 2 . This is of opposite sign compared 
to the values obtained by Oh et al£^ } and also larger 
than their values of 72 and 130 Oe/V 2 . On the other 
hand, in the P configuration the resistance of the device 
is 1.7 times lower than in the AP configuration, so for 
the same applied bias voltage, the power dissipated by 
Joule heating is a factor of 2.8 higher in the P config- 
uration than in the AP one. The Joule heating is not 
enough to explain the softening of a FL mode due to re- 
duced magnetization density, which would in turn reduce 
the restoring torque provided by the demagnetizing fields 
(N y — N x )Ms- The Curie temperature of the FL is close 
to 1000 K, and the temperature necessary to significantly 
reduce the magnetization density would have destroyed 
the tunnel junction. However, if the FLP1 is actually 
a coupled mode in which both FL and RL oscillate, but 
with different amplitudes, this can both explain the lower 
peak power spectral density of this mode as well as the 
bias dependence of the mode. Any in-phase correlated 
motion between the FL and the RL will necessarily re- 
duce the resistance fluctuations of the mode, since they 
are proportional to the cosine of the angle between the 
RL and FL magnetization directions, and therefore also 
the power spectral density of the mode. A main restoring 
torque of such a coupled RL-FL mode will be provided 
by the exchange bias, and a reduction in exchange bias 
due to Joule heating will soften the mode and can possi- 
bly lead to the observed (weak) parabolic dependence of 
this mode's frequency. 

V. SUMMARY AND CONCLUSIONS 

In summary, we have measured the resonance fre- 
quency of a nominally circular MgO tunnel junction with 



approximately symmetric electrodes, both as a function 
of applied external magnetic field and as a function of 
bias voltage, and analyzed the possible dependence of 
the perpendicular spin torque effective field as a function 
of bias voltage for different, fixed magnetic fields. At low 
magnitudes of magnetic field and low bias voltage, the 
dominant low- frequency mode (FLAP1) is a FL bulk- 
like mode. At larger fields in the AP orientation, there 
is a crossing of this FL branch and a RL branch, and the 
low-frequency mode is a bulk-like RL mode (RLAP1). 
At negative fields in the P orientation, the mode appears 
to be a FL mode (FLP1), that could be more localized 
to an edge than the FLAP1 and RLAP1 modes. The 
frequency of the FLAP1 mode at fixed magnetic field 
has a bias voltage dependence that is linear in bias volt- 
age, which leads us to conclude that the perpendicular 
spin torque effective field is linear in bias voltage with 
a value of about 37 Oe/V, or 0.55 x 10~ 6 Oc(cm) 2 /A. 
This is in qualitative agreement with previously observed 
linear dependence of the perpendicular spin torque ; 18 i 23 
although the magnitude of the effect differs, but in dis- 
agreement with the quadratic dependence observed by 
other groups i 12 ' 19 i 20 ' 24 The quadratic behavior of the fre- 
quency of the RLAP1 mode that we observed is entirely 
due to Joule heating, which reduces the exchange bias 
and consequently the restoring torque for this mode. Any 
(weak) linear or quadratic dependence on the perpen- 
dicular spin torque effective field cannot be extracted 
beyond the experimental uncertainties. For the FLP1 
mode, a much weaker quadratic dependence of the reso- 
nant frequency on bias voltage was observed. We do not 
believe this quadratic dependence is due to a quadratic 
dependence of the perpendicular spin torque effective 
field upon bias voltage. The reason for this is that the 
sign is opposite to that obtained by Oh et aS^ and the 
effective field would be about a factor of three larger in 
magnitude. We speculate that this bias dependence of 
the mode frequency is due to a weak dependence on ex- 
change bias of what is actually a (weakly) coupled FL-RL 
mode. 
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